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Abstract: We consider an atom inside a ring cavity, where a plane-wave cavity field together 
with an external coherent laser beam induces a two-photon Raman transition between two 
hyperfine ground states of the atom. This cavity-assisted Raman transition induces effective 
coupling between atom’s internal degrees of freedom and its center-of-mass motion. In the 
meantime, atomic dynamics exerts a back-action to cavity photons. We investigate the 
properties of this system by adopting a mean-field and a full quantum approach, and show 
that the interplay between the atomic dynamics and the cavity field gives rise to intriguing 
nonlinear phenomena. 
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1. Introduction 

When Jaynes and Cummings first studied the time evolution of a two-level atom in an electromagnetic 
field in a fully quantized way in 1960s [1], experimental realization of this ideal theoretical model was 
out of reach. It was made possible only with the advent of one-atom masers in late 1980s, by Rempe, 
Walther and Klein [2], who experimentally studied the interaction of a single atom and a single resonant 
mode of electromagnetic field in a cavity. Jaynes-Cummings model (or J-C Model) serves to bridge 
our understanding of the relationship between quantum theory of radiation and semi-classical theory 
of atom-light interaction, and has become one of the most important models in quantum optics and 
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cavity electrodynamics (CQED). In the experiment of Referenee [2], a beam of atoms was used and the 
atom-light interaetion was studied during the transient time when the atoms pass through the eavity. The 
intensity of the atomie beam is suffieiently low sueh that at any given time, no more than one atom is 
inside the eavity. The field of CQED was further advaneed by putting a single trapped eold atom [3], 
and more reeently a eondensate of ultraeold atoms [4-6,8], inside an optieal eavity. In an optieal ring 
eavity, through eoherent eontrolling the dark and bright states, the vaeuum Rabi splitting and optieal 
bistability of eavity multi-wave-mixing proeess has been studied [9]. In addition, bright-state polaritons 
of four-wave mixing and six-wave mixing signals ean be parametrieally amplified [10]. In the J-C model, 
and other related CQED systems, the foeus is the interaetion and mutual influenee between the eavity 
mode and the atomie internal degrees of freedom. The external degrees of freedom of the atom, i.e., its 
eenter-of-mass (COM) motion, is generally negleeted. Unlike “hot” atoms, however, eold atoms’ COM 
motion in general ean no longer be negleeted in this “atom -i- eavity” system, as the COM momentum 
of a eold atom will be signifieantly affeeted by photon reeoil from emission and absorption of even a 
single photon. Therefore in a more eomplete deseription of the eavity system, one needs to take into 
aeeount the interplay among the eavity photons, and both the internal and external atomie degrees of 
freedom. Eurthermore, when more than one atom is inside the eavity, one should also aeeount for the 
eavity photon mediated long-range interaetion between atoms. All of these make the “eold atoms -i- 
eavity” system extremely rieh and interesting, and truly represent a new frontier in both CQED and eold 
atom researeh. 

As the same photon affeets both the internal states (via indueing a transition between different states of 
the atom) and the external COM motion (via photon reeoil) of the atom, it naturally induees a eoupling 
between the two atomie degrees of freedom. Sueh spin-orbit eoupling (SOC) in eold atoms has been 
realized in both bosonie [11,12] and fermionie systems [13,14], and has attraeted tremendous attention 
in reeent years [15]. In praetiee, to avoid spontaneous emission, SOC is indueed between two hyperfine 
ground states of an atom via a pair of Raman laser beams. Due to its non-Abelian nature, SOC not only 
signifieantly affeets the physies of a single atom, but, perhaps more importantly, also profoundly ehanges 
the properties of a many-body system. It is an essential ingredient underlying sueh diverse phenomena 
as topologieal supereonduetors/insulators, Majorana and Weyl fermions, spin-Hall effeets, etc. [16-20]. 

So far, all the experimental realization of photon-indueed SOC in eold atoms employs elassieal laser 
fields to generate the Raman transition. Here the parameters of SOC are determined onee the laser 
parameters (e.g., intensity, frequeney, etc.) are fixed. In a reeent work [21], we propose to replaee one 
of the Raman beams by a plane-wave eavity mode supported by a ring eavity. Around the same time, 
Mivehvar and Eeder eonsidered a situation where both Raman beams are provided by a eavity ([22]). 
In this seheme, as the baek-aetion from the atoms affeets the eavity photons, the SOC thus indueed 
beeomes nonlinear and dynamieal. By employing a simple mean-field approaeh, we demonstrated [21] 
that sueh a system indeed possesses interesting nonlinear properties even with just one single atom inside 
the eavity. Eor example, the eavity-assisted SOC dramatieally modifies the atomie dispersion relation, in 
partieular, with the emergenee of loop struetures under eertain eireumstanees. Several other groups have 
reeently investigated the eavity-assisted SOC in many-body atomie systems [23-26]. 

In the present work, we first briefly review our previous proposal [21], and then theoretieally explore 
the full quantum meehanieal treatment beyond the mean-field formalism, and finally investigate the 
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correspondence between the quantum and the mean-field treatment. The quantum treatment is carried 
out by solving the Master equation for the total density operator, from which we can derive various 
quantities of interest, i.e., the cavity photon statistics, the degree of entanglement between the atom and 
cavity field, etc. These two different approaches provide a deeper understanding to this intriguing system. 

The article is organized as the following: After briefly reviewing key ideas of our previous work 
and the mean-field approach in Section 2, we develop the full quantum mechanical formalism to the 
physical system of interest in Section 3 and discuss about the intimate correspondence between the two 
in Section 4, and finally conclude in Section 5. 

2. Model Setup and Mean-Field Formalism 

As shown schematically in Figure 1, we consider a single atom being confined by a single-mode 
unidirectional ring cavity, whose cavity mode together with an additional coherent laser beam form a 
pair of Raman beams that induces transition between two hyperfine atomic ground states denoted as 
I t) and I J,), while transferring recoil momentum ±2hqrZ to the atom. The ring cavity has a resonant 
frequency of cUc, decay rate k, and is pumped by an external laser field with frequency cUp and pumping 
rate Cp. The Hamiltonian under the rotating wave approximation can be written as [21], 

'HeS ^ ^ 

a=t,; 

+ itp{c^ - c) - (1) 

where, for simplicity, we only consider the atomic COM motion along the ^-axis, which is the direction 
of the photon recoil. In Hamiltonian Equation (1), k is the atomic COM quasi-momentum (we take 
^ = 1) along the z-axis, c is the cavity annihilation operator, r[)(j(z) (cr =1, ],) is the atomic operator, 
= ±1, 25 represents the two-photon Raman detuning, 5c = cUp — cUc is the cavity-pump detuning, 
and O denotes the atom-cavity coupling strength (or more specifically, the Raman coupling strength per 
cavity photon). Finally we will treat the cavity decay phenomenologically, which amounts to adding a 
non-Hermitian term —inc^c in the above effective Hamiltonian. Note that for a quasi-momentum k, the 
real momentum for an atom is spin-dependent: it is A: -f qr for | ]') state and k — qr for | ],) state. 

From Hamiltonian (1), one can easily obtain the Heisenberg equations of motion for both atomic and 
cavity fields. In this work, we only consider a single atom inside the cavity. Hence we rewrite the atomic 
field operators as first quantization wave functions —)■ 4’( 7 ( 2 ^)- Assuming spatial homogeneity, we 

further take the plane-wave wave function for the atomic modes 4 >(r( 2 :) = with the normalization 

condition |(p^p -f |(p^p = 1. To proceed further, we adopt the mean-field approximation by replacing 
the photon field operators by its respective expectation value: c —)■ c = (c). The validity of this 
mean-field approximation lies on two assumptions: (1) the cavity field can be described as a coherent 
state; (2) the atomic field and the cavity field have negligible entanglement. We shall come back to 
these two assumptions later when we compare the mean-field results with results obtained from the 
beyond-mean-field quantum Master equation. 




-f 2lXaqrk 
2m 


+ oc^S r[)o-(2;) 




^iz)^l{z)c + h.c. 
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Figure 1. Schematic diagram of the cavity-assisted spin-orbit coupled system. Inside the red 
dashed box, we show the level diagram of the atom and the light field configuration where 
cUp represents the cavity field and ujr is an external coherent laser beam. Zp and k are the 
cavity pumping rate and decay rate, respectively. 


Within the mean-field formalism, the steady-state solution for the photon field is obtained by taking 
the time derivative of the photon field to be zero, from which we obtain: 


c = 


Zp-i^(pl(p^ 


K — iSr 


( 2 ) 


Inserting this into the equations for atomic fields, we obtain the coupled nonlinear time-dependent 
equations for the two spin components as. 


I , c\ ^eff 


(3) 

(4) 


where fleff = flc = effective Raman coupling strength between two atomic states. 

The fact that fles depends on the atomic field itself is a manifestation of the non-linearity arising from 
the atomic back-action to the cavity field. 

For a given atomic quasi-momentum k, we define eigenstate and eigenenergy as the solution of the 
time-independent version of Equations (3) and (4), by replacing i{d/dt) with e{k). After some lengthy 
but straightforward algebra, we find that e{k) obeys a quartic equation in the form of 


4e^ + Be^ + Ce^ + De + E = 0 


(5) 


The derivation of the above equation and the explicit expressions of the fc-dependent coefficients B, C, 
D and E can be found in our previous work [21]. 

We can gain some insights about the general structure of the dispersion relation e{k), e.g., the 
degeneracy condition and the appearance and disappearance of the loop. In the case of vanishing 
two-photon detuning (i.e., 6 = 0), simple analysis shows that there should be a total of four regimes. 
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as shown in Figure 2a. In region I, the two dispersion branehes are gapped, and the lower braneh has a 
double degenerate minima, as shown in Figure 2(bl). This dispersion eurve strueture is very similar to 
the ease when both Raman beams are provided by elassieal eoherent laser fields (we shall refer to this as 
the “elassieal ease”) and the Raman eoupling strength is small. In region II, as shown in Figure 2(b2), 
the two dispersion branehes are still gapped, but the lower braneh has a single minimum. This is similar 
to the elassieal ease with a large Raman coupling strength. Regions III and IV do not have analogs in 
the classical case. Region III features a loop structure, as shown in Figure 2(b3), whereas in Region IV, 
the loop dissolves but the two dispersion branches becomes gapless at A; = 0, as shown in Figure 2(b4). 
In the looped region, the quartic Equation (5) yields four real roots. It can be shown [21] that this requires 
the coupling strength to satisfy 


=^tp<^< 4£pVl + (6e/K)2 = 


The two dashed lines in Figure 2a represent the two critical values and respectively. 





k/qr 


Figure 2. Single particle eigen-energy spectrum “phase diagram”. The dispersion curve 
is generally catagorized by four regions, represented from I to IV in (a). From (bl) to 
(b4), we fix £p = k. In region I, the dispersion has double minima as shown in (bl) with 
= 0.03 k; region II is enclosed by the red solid curve in (a) and we show the typical point 
in (b2) (f2 = n) where only single minimum exisits in the lower helicity branch; region III 
is enclosed by the black dashed lines in (a) and it is a region where loop structure emerges, 
as in (b3) with = 5 k; finally, in region IV we recover the double minimum dispersion 
although it’s different from region I by closing the gap at /c = 0, as in (b4) with Vt = S n. 
Throughout the paper we fix 6c = k and 6 = 0, and adopt a dimensionless unit system with 
Ti = m = K = l. A typical value for k is 27r x 1 MHz, and we choose qr = 0.22 in our 
dimensionless units (based on a realistic experimental parameter estimate). 


It is instructive to examine how the effective Raman coupling strength fleff behaves as a function of 
In Figure 3a, we plot |f2eff| for the lowest dispersion branch as a function of for different k values. 
Note that |f2effP = ^^’^photon where riphoton = |cp is the steady-state cavity photon number. A few 
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remarks can be made based on this plot. First, the fact that fleff is different for different k clearly shows 
the influence of the atomic COM motion on both the atomic internal dynamics and the cavity photon 
number. Second, |Oeff| is not a monotonous function of O. For given k, |Oeff I increases with O linearly 
for small O. This can be intuitively understood as follows. At such weak atom-photon coupling, the 
back-action from the atom to the cavity photon is negligible. The number of cavity photons riphoton is 


roughly given by riphoton ^ Uq = 
atom is absent. As a result, we have |f2eff 


K, — iSr 




, where uq is the number of cavity photons when the 
kly/rio which is independent of the atomic quasi-momentum 
k. On the other limit, when O is very large, the strong atom-cavity coupling strength significantly detunes 
the cavity away from resonance and the cavity photon number riphoton, and hence |Oeff|, decreases as a 
function of O. Such a non-monotonous behavior of Oeff is a unique feature of the cavity system and a 
direct manifestation of the non-linearity of the system arising from the back-action of the atom on the 
cavity photon. 



0/b: 



Figure 3. (a) Effective Raman coupling |f2eff| is plotted as a function of atom-photon 
coupling strength kl for different k values, 0, qr, 10 qr for blue dash-dot, red solid and black 
dashed lines. We observe that |f2eff| does not monotonically increases with kl but rather 
peaks at an intermediate value, then approaches zero in the large limit. Figure (b) shows a 
comparison between critical boundary of region I and II (red solid curve) and the analytical 
result (blue dashed line) given in Equation (6). At large tp limit, the two results match 
asymptotically well. 


From the above analysis, it should also become clear that when the effect of the back-action is weak 
(which occurs when is small and/or Zp is large), we should recover the properties of the classical case. 
In particular, in the classical case, the lower dispersion branch change from two degenerate minima to 
a single minimum when the Raman coupling strength exceeds a critical value. Using our notation, this 
occurs when |Uefr| exceeds the critical value AEr where Er = ql/{2m) is the photon recoil energy. 


For weak atom-cavity coupling, |f2eff I = = U 




. Hence the critical value of U is given by 


U = 4U 




( 6 ) 
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In Figure 3b, we plot this eritieal value (blue dashed line) as a funetion of eavity pump rate £p and 
eompare it with the numerieally determined lower boundary (red solid line) between region I and II of 
Figure 2a. The two eurves overlap with eaeh other when £p inereases. Therefore, as we have expeeted, 
in the limit of weak atom-eavity eoupling and strong eavity pumping, we fully reeover the elassieal ease 
where the SOC is indueed by two elassieal laser beams. 

3. Master Equation Approach: Full Quantum Mechanical Treatment 

The above diseussion is based on the mean-field approaeh where the eavity field is replaeed by 
a c-number that represents the photon amplitude. This mean-field treatment relies on two implieit 
assumptions: (1) the atom-photon eorrelation is negligible, and (2) the photon field ean be well 
approximated by a eoherent state. In order to examine the validity of these assumption, and henee 
the validity of the mean-field approximation, we now turn to a full quantum treatment based on the 
Master equation: 

p=^[H.n,p]+C[p] (7) 

Here p is the total density operator of the eoupled atom-eavity system, the effeetive Hamiltonian "Heff is 
the same as in Equation (1). The dissipation arising from eavity deeay is modeled by the Liouvillean 
term in the standard form of Lindblad super-operator [27,28], 

C[p] = K{2cpc^ — c^cp — pc^c) (8) 

Again, due to spatial homogeneity, we deeouple momentum eigenstates by taking the plane-wave ansatz 
for the atomie modes as 'ii>a(2:) = As there is no eoupling between atomie operators with 

different k, we ean work in the subspaee for a fixed value of k. Here we explieitly write the eommutator 
in Equation (7), for a given atomie quasi-momentum k, as, 

[w=,(fc). p] = (^ + ^+s) ('h'i’rf - I'i’t) + (Is ■ ^ ■ ''’*’1**) 

+itp (c^p — cp — pc^ -f pc) — 6c (c^cp — pc^c) (9) 

Note that if the photon reeoil = 0, whieh oeeurs when the eavity mode and the external laser beams 
are eo-propagating, the COM kinetie energy terms k"^/{2m) ean be gauged away after a simple gauge 
transformation. Our model is then redueed to the J-C model and the atomie COM motion does not play 
a role. To solve the Master Equation (7), we ehoose our basis states as direet produet states of photon 
Eoek state In) and atomie internal state |cr =t, i): |ri; o) = |n) 0 |cr), where non-negative integer n 
denotes eavity photon number. Our goal is to ealeulate the entire matrix elements of the density operator 
under this set of basis states, denoted by (m; (y\p\n] o') = p////. We found the governing equation for the 
matrix element ean be written as. 
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dt 


\2m 


m 


(ja 

mn 


n 


+ IPm+ln + ^ols/^Pm-ln ^tr'^Vn + lp^^_|_^ ^0'4V^Pmn-l) 

+ Ep (^V^Pm°-ln ~ Vm + lPm‘+ln + \AiPmn-l ~ \/n + 

+ z6c (m - n) + K (^27 m + lVn + lp^+i„+i - (m + 


( 10 ) 


where ct represents the flip-spin value, i.e., f =1 and | =t. 

With Equation (10), we ean study the dynamieal evolution of the density operator p for a given 
initial state. Obviously, we need to introduce a sufficiently large photon number cutoff. Once we obtain 
the density operator, all relevant physical quantities can be readily calculated. An example is given 
in Figure 4, where we show the time evolution of the cavity photon number n = Tr[ph] = Tr[pc^c] 
for the initial state |0;t)- The three different curves in Figure 4 correspond to different atomic 
quasi-momentum k. 



tK 


Figure 4. Time evolution of cavity photon number. The initial state is given by |0;t) and 
we consider the same parameter as in Figure 5a with k/qr = —10, 0, 10. The steady-state 
values, obtained in the long time limit as shown here, correspond to red dashed lines at 
corresponding k values in Figure 5a. 


As evidenced in Figure 4, due to the presence of cavity decay, a steady state will eventually be reached. 
Fet us now focus on the properties of the steady state. The steady-state density operator matrix elements 
can be obtained by equating the RHS of Equation (10) to zero. The red dashed lines in Figure 5a-c 
represent the steady-state photon number as functions of atomic quasi-momentum k. The horizontal 
arrows in the plot represent the cavity photon number by setting = 0, in which case our model 
reduces to the J-C model and all physical quantities become /c-independent. 
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Figure 5. (a)~(c) Photon number obtained from the mean-field approach (solid curves) and 
from quantum mechanical Master equation approach (red dashed curves). From (a) to (c), 
= 3 K, 5.6 K, 6 K. The color on the solid curves represents the normalized decay rate 7/f2 
of unstable mean-field states. The black color represents stable mean-field states. We have 
used tp = K, and other parameters are the same as in Figure 2. (d)~(f) Corresponding 
photon number fluctuation (blue solid curve) and negativity (green dashed line) obtained 
from the quantum approach. The parameters used here are the same in (a)~(c), respectively. 
The horizontal arrows indicate results from the J-C model by taking = 0. 


To have a better understanding of the photon statistics, we study the steady-state photon number 
fluctuation. Specifically, we calculate the normalized photon number fluctuation defined as 

((An)^) (h^) - (h)^ 

{n) (n) 

where the expectation values of the operators are obtained with the help of the steady-state density 
operator. For a coherent cavity field, the photon fluctuation is Poissonian and we have A„ = 1. The solid 
curves in Figure 5d-f represent A„ (left vertical axis) as functions of k, and the horizontal arrows 
pointing to left give the values of A„ from the J-C model by setting qr = 0. For the parameters we 
have used, we note that the J-C model always predicts a super-Poissonian photon statistics, whereas our 
model gives super-Poissonian photon statistics only for small atomic quasi-momentum, but Poissonian 
statistics as k/qr —)■ ±oo. 

Last but not least, to characterize the correlation between the atom and the cavity field, we investigate 
the so-called negativity [29] which measures the degree of entanglement for a mixed state system. 
The negativity is defined as Afip) = (||p^'^||i — l)/2, where is the partial transpose of the density 
operator with respect to either the atom subsystem or the cavity subsystem, and | | |i denotes its trace 

norm with the definition | |A| |i = Tr[\/ A^A]. A negativity of zero indicates that the two subsystems (the 
atom and the cavity, in our case) are not entangled, whereas a positive negativity means that finite degree 
of entanglement is present. The dashed curves in Figure 5d-f represent the negativity (right vertical axis) 
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in the steady state as funetions of k, and the horizontal arrows pointing to right give the values of the 
negativity from the J-C model by setting = 0. One ean observe that for the ehosen parameters, the J-C 
model always prediets a finite degree of entanglement between the atom and the eavity field. By eontrast, 
the degree of entanglement in our model weakens when k/qr —)■ ±cx). 

4. Discussions 

In the previous two seetions, we have presented two different methods for studying the system. We are 
now in a position to diseuss their eonneetions. 

In Figure 5a-e, in addition to the steady-state photon number obtained from the quantum treatment 
(dashed eurves), we also plot the photon number riphoton = |cp obtained from the mean-field approaeh 
(solid curves), with c given in Equation (2). In the quantum treatment, the steady-state density matrix is 
obtained by solving a set of linear equations. For a given k, the solution is unique. Hence we only get 
one steady-state photon number for a given atomic quasi-momentum. On the other hand, the mean-field 
treatment allows multiple steady-state solutions corresponding to different real roots of the quartic 
equation Equation (5). Hence a single k value is associated with more than one steady-state photon 
number. However, due to the non-linearity intrinsic in the mean-field method, not all steady-states are 
dynamically stable. A straightforward stability analysis allows us to quantify the dynamical stability of 
the mean-field states, as we did in Reference [21]. The stability information of the mean-field states 
are encoded as the color value in the solid curves. Those stable states are represented by black color, 
while any color other than black indicates an unstable state, and the color value represents the decay 
rate (see the colorbar) of the corresponding state. From Figure 5a-c, we clearly see that at large atomic 
quasi-momentum k/qr ^ ±oo, the Master equation result overlaps with the stable mean-field branch; 
while at small \k\, the quantum result deviates significantly away from the mean-field solution. 

The agreement for large \k\ and the discrepancy at small \k\ are both consistent with the results of 
the negativity and photon number fluctuations as presented in Figure 5d-f: At large \k\, the negativity 
is small {i.e., atom-cavity entanglement is weak) and the photon number fluctuation tends to Poissonian 
{i.e ., the photon field is well approximated by a coherent state), this is exactly the regime where we expect 
the mean-field approximation is valid. By contrast, for small \k\, the quantum calculation indicates that 
there is non-negligible entanglement between the atom and the cavity field, and the cavity field itself 
cannot be assumed as a coherent state. Hence the mean-field assumption is no longer valid. 

The reason why mean-field approximation only works for large | A:| is actually rather simple. Consider 
a Raman transition process where the atom jumps from | j') to | ],). The quasi-momentum k does not 
change during this process, however the real momentum changes from k + qr to k — qr. Therefore 
the effective two-photon Raman detuning is not just 26, but 26 -f 2qrk/m, where the additional term 
comes from the difference of the kinetic energies for different pseudo-spin state |cr). In other words, the 
SOC renders the two-photon detuning momentum-dependent. In the examples we presented in this 
work, we have taken 6 = 0. Hence the Raman transition is only near-resonant for small \k\, and 
becomes far off-resonant for large \k\. Therefore, for large \k\, the atom-photon coupling, and hence 
the atomic back-action to cavity, are weak. This explains why the mean-field assumption becomes valid 
in this regime. 


Atoms 2015, 3 


192 


5. Conclusions 

In this work, we have studied spin-orbit coupled cold atoms inside a ring cavity system, employing 
both the mean-field theory and the full quantum mechanical Master equation approach. By treating 
both light and atom on equal footing and seeking the self-consistent solution in both approaches, we 
have found: (1) cavity-assisted SOC dramatically modifies atomic dispersion relation, (2) intriguing 
dynamical instabilities exist in the system, (3) atom’s back-action onto cavity field also leads 
to non-trivial atom-photon coupling that is fundamentally different from either the system with 
classical-laser induced SOC in the absence of the cavity or the J-C model where the atomic COM motion 
is neglected (Le.,hy taking = 0). We have also explored correspondence and discussed the connection 
between the mean-field and the quantum approaches. The two distinctively different approaches provide 
us with a deeper understanding and complementary insights into this system. We conclude that the 
synthesis of cavity QED and SOC is not a trivial combination and interesting new physics emerges in 
this setting. 
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